Geometric effects on T-breaking in p + ip and d + id superconductors 
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Superconducting order parameters that change phase around the Fermi surface modify Josephson 
tunneling behavior, as in the phase-sensitive measurements that confirmed d order in the cuprates. 
This paper studies Josephson coupling when the individual grains break time-reversal symmetry; the 
specific cases considered are p±ip and d±id, which may appear in S^RuCU and Na a; Co02-(H20)j / 
respectively. T-breaking order parameters lead to frustrating phases when not all grains have the 
same sign of time-reversal symmetry breaking, and the effects of these frustrating phases depend 
sensitively on geometry for 2D arrays of coupled grains. These systems can show perfect supercon- 
ducting order with or without macroscopic T-breaking. The honeycomb lattice of superconducting 
grains has a superconducting phase with no spontaneous breaking of T but instead power-law corre- 
lations. The superconducting transition in this case is driven by binding of fractional vortices, and 
the zero-temperature criticality realizes a generalization of Baxter's three-color model. 

PACS numbers: 74.50+r 74.72-h 



I. INTRODUCTION 

Several materials are now believed to support su- 
perconducting states where the order parameter A(k) 
is complex-valued and hence breaks time-reversal sym- 
metry. This class includes the recently reported 
superconductor^ Na^CoG^^EbO),,, which is a doped 
Mott insulator since C0O2 is a layered triangular- 
lattice spin-half antiferromagnet. RVB-type mean-field 
theories^*^ predict that the pairing symmetry should be 
d±id, and successful observation of this state would lend 
strong support to the RVB picture of unconventional su- 
perconductivity. T-breaking superconducting pairing is 
also believed to occur in the layered material S^RuO^ 5 , 
where triplet Cooper pairs have zero spin along the per- 
pendicular direction (S z = 0) and two degenerate orbital 
states p±ip with angular momentum ±h in the perpen- 
dicular direction. 

A major experimental success in the study of the 
cuprates was the use of phase-sensitive measurement to 
confirm the d x 2_ y 2 pairing symmetry^* 7 ". In these exper- 
iments the coupling between the superconducting order 
parameter and the underlying crystalline lattice is used 
to design a frustrated loop of weak links. (Each frus- 
trated loop is associated with a "flux" which is the phase 
of the product of the Josephson coupling constant around 
the loop.) Similar experiments may constitute a powerful 
test of p + ip and d + id symmetries (section [nj. 

The main purpose of this paper is to demonstrate the 
rich physics of arrays of Josephson junctions made up 
of superconductors with T-breaking order parameters. 
Throughout this paper we shall focus on two dimensional 
arrays. In particular we shall show that even in the ab- 
sence of external magnetic field there is a effective flux 
associated with each elementary loop that frustrate the 
Josephson energy. This flux is determined by the Ising- 
like time-reversal order parameters in the grains along 
that loop. The Ising order parameter specifies which of 



two degenerate states related by time reversal exists in a 
particular grain: for example, d + id (a = +1) or d — id 
((j = —1). The classical Hamiltonian of such a system 
takes the Higgs-gauge form 



H 



HG 



= -Ej 



<ij> 



cos 



<j)j - ^-(er^CTj))- 



(1) 



In the above 4n is the phase of the superconducting order 
parameter in the ith grain and Oi is the Ising-like time 
reversal order parameter in that grain. The gauge field 
Aij is given by 



As 



has 



(2) 



where k = 1 for p ± ip and k — 2 for d ± id, and 0i is 
the angle between the crystalline y-axis of the ith grain 
and the normal of the interface from i to j. The ex- 
istence of these discrete (Ising) variables as well as the 
continuous superconducting phase variables, and the fact 
the former affects the latter as a gauge field, is the main 
difference between the problem we study here and the 
usual frustrated/unfrustratcd Josephson junction arrays. 
Frustrated s Josephson arrays have recently attracted at- 
tention as a model for quantum computation 8,9 . Direct 
coupling between local impurity spin variables and a su- 
perconducting order parameter is discussed iniS. 

Unlike in the real d or s case, Ginzburg-Landau the- 
ory for T-breaking superconductors^ does not correctly 
give even the phase diagram of Josephson arrays. For 
example, the coupling between superconducting phase 
and array geometry can induce classical discrete Z2 or 
Z3 topological order, where a fractional vortex of the 
superconducting phase is bound to a defect of the local 
Ising variables that carries fractional Ising quantum num- 
ber. This is analogous to the Quantum Z2 topological 
order that has attracted considerable attention recently 
because of its appearance in quantum dimer and other 
modela 12 ' 13 . 
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Many of the nontrivial properties of Josephson arrays 
of T-breaking superconductors are connected to models 
of frustrated magnets. For example, the Ising configu- 
rations {<7i} that leave the interaction among the phase 
variables unfrustrated are equivalent to the ground-state 
spin configurations of certain frustrated magnets. On 
the honeycomb lattice, the constraint of no frustration 
for p ± ip or d ± id superconductors becomes the same 
as that in the "three-coloring problem"—, studied earlier 
because of its relation to classical frustrated antiferro- 
magnets on the kagome latticei^*iSiiLi2iiS. As in the XY 
antifcrromagnet on the kagome, the zero-temperature su- 
perconducting phase has long-range order of ip 3 but not 
ip, where ip is the Cooper pair operator. That is, the 
phase of A(ko) for some fixed ko is not uniform through- 
out the array in the low-temperature phase, but varies 
only by integer multiples of 2tt/3. 

Another example of our results is that, on the square 
lattice, the Ising configurations that do not cause phase 
frustration for p + ip superconductors are the same as 
the ground state gauge field configurations of a Z 2 gauge 
theory. The low-temperature state is like a nematic in 
that ip 2 but not ip has an expectation value. 

This paper is organized as follows. In section II we cal- 
culate the Josephson tunneling amplitude between two 
T-breaking superconductors and from microscopic con- 
siderations. In sections III we discuss some general fea- 
tures of the Higgs-gauge model given in Eq. J[J. In sec- 
tions IV through VI we study the novel physics of Eq. (Q) 
on different lattices. This corresponds to Josephson ar- 
rays of T-breaking superconductors in different geome- 
tries. Section VII discusses the entropic contribution to 
the interaction free energy between fractional vortices. 
Finally in section VIII we summarize the main conclu- 
sions and some connections to other recent work. 



II. THE JOSEPHSON COUPLING BETWEEN 
TWO COMPLEX SUPERCONDUCTORS 



Consider the tunneling Hamiltonian between two su- 
perconducting domains 1 and 2 



H t = - (*( fc lh /c i' /c l) ct i;fc ll ,fc_ L c 2;fc||,fc^ +h.c.j (3) 

fell ,k± ,k' ± 

Here the momentum parallel to the interface (fc|| ) is con- 
served by the tunneling process, while the perpendicular 
momentum (k±) is not, as appropriate for a sharp bound- 
ary. 

Let us assume the crystalline y-axis of the two domains 
are oriented at angles 0i and 02 with respect to the inter- 
face normal. The Josephson coupling between these two 
domains can be calculated from second-order perturba- 



tion theory. The result is 

£i2(0i,0 2 ; tufa) =- 



E 



fen ,k± 



|*(fcii,fc±,fcl) 

2A 



A 1 (A||,*x)AS(fc||,tf L ) + && ,(4) 



where A is the superconducting order parameter. 
In deriving the above we have used the fact that 
t(-k\\,-kx,-k' ± ) = t*(k\\,k±,k' ± ). 

Now we specialize to the d ± id case. With respect 
to the individual crystalline axis (x, y) the d + id order 
parameter has the form 



A(k x ,k y ) = /\ e^{k x +i(jk v Y 



(5) 



where <f> is the phase of the superconductor order param- 
eter and a is the time reversal order parameter 2 ^. When 
transformed to the interface coordinate system we have 



Ai(fc||,fc| 
A a (fc||,fcx) = A . 



A n e 2i<Tiei < 



1,1 (k 



,2i<j 2 2 i<t>2 



(fell 



i<Jik^)' 2 
ia 2 k±) 2 . 



Substituting Eq. © into Eq. (@J we obtain 

E12 = -Ej cos(<f)i - <f> 2 + ^12(01,0-2)) 

Al2(<7l, (T2) = 2(7101 — 2<7 2 02. 

In the above 

|t(fc||,fcx,fcx)| 2 



(6) 
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E, 



- A 2 
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E 

fell ,k±,k' 



2A 



-(fc|| + ia]k\_y 



x(fc|| - ia 2 k' A _f 



(8) 



Note that when \t(k\t , k±, k' ± )\ 2 is even upon reversing k± 
and k' ± (as should be the case for local tunnelling) Ej is 
independent of a\ and a 2 ■ 

Similar expressions can be obtained for p + ip pair- 
ing except in this case Ai 2 (ai,a 2 ) = cr 1 1 — C202- For a 
straight boundary, 9 2 = 0i+7r. This ir phase between op- 
posite orientations is unimportant for even angular mo- 
mentum order parameters (e.g., d+id) but significant for 
odd angular momentum order parameters (e.g., p + ip). 



III. JOSEPHSON JUNCTION ARRAYS OF 
T-BREAKING SUPERCONDUCTORS 

We now show from Eq. (JJJ that the classical Hamilto- 
nian for an array of p + ip or d + id Josephson junctions 
is given by the gauge theory Eq. Q and Eq. (J2J. The 
Josephson Hamiltonian in terms of gauge- invariant phase 
differences across each junction is 



H = —Ej y^cosc 

(ij) 



(9) 



These phase differences are related to the superconduct 
ing phases on the grains by 4>ij — <Pi ~ 4>j 



Aij with 



Aij = koidi — kcfjOj. The partition function of the sys- 
tem is 



* = II£ 



. e P E -jT.ij cos(<p i - cjij -Aij (a i, a j)) ^ ^q-j 



This theory has a classical compact U(V) gauge invari- 
ance under 



Xr ~ Xj- 



(11) 



Standard duality transformational* 2 ^ 3 " applied to the 
above model gives the vortex partition function 

{&j}{n a } 

In the above a, b label sites of the dual lattice, the n a are 
integers, and the flux is 



(13) 



<ij >€a 



A 2 

The fugacity of vortices is e~ cUa and G a b is a loga- 
rithmic interaction between the vortices. In this form, 
the problem becomes a 2D vortex gas with annealed flux 
variables. Each flux configuration is weighted by a purely 
cntropic term that measures how many Ising configura- 
tions correspond to it. This type of problem has been 
encountered previously in Ref. 23 . 



IV. TRIANGULAR LATTICE: CONFINEMENT 

As a first step in understanding how the gauge phases 
in Eq. affect macroscopic ordering, we consider a reg- 
ular triangular lattice of d + id or p + ip grains (Fig. 
For an elementary plaquette with three corners 1, 2, 3, it 
is simple to show that the sum of A+j gives rise to flux 



k 

3 (ai 



cr 2 + cr 3 ) 



mod 1. 



(14) 



The constant k is 1 for p+ip and 2 for d+id. The ground 
state configuration of {ai} satisfies the constraint that <i> Q 
be an integer for all plaquettes (a runs through a honey- 
comb lattice). The only way to satisfy this constraint is 
for all <7j to take the same value +1 or —1. 

Flipping a single a creates six neighboring plaquettes 
with fractional $ a = ±2fc/3, so the total flux is inte- 
gral. It is possible to create local fractional fluxes at 
domain wall "kinks" : these are corners in a line sepa- 
rating a region of Ising spin +1 from one of Ising spin -1 
(Fig. 121) • These fractional fluxes are connected by a defect 
line which costs finite energy per unit length. The cn- 
tropic term is quite strict for the triangular lattice, since 
any flux configuration corresponds to at most two spin 
configurations. As a result, the fractional vortices are 
confined. To summarize, the triangular case has perfect 




FIG. 1: A triangular lattice of d ± id superconducting grains 
connected by Josephson tunneling. The state of grain j is de- 
termined by an Ising variable Oj = ±1 and a superconducting 
phase. 




FIG. 2: Domain wall kink on the triangular lattice: there is 
a localized fractional flux <l?o/3 from the plaquette inside the 
dotted triangle. Here solid (open) circles correspond to spin 
up (down). 



order of both the phase and Ising variables at T = O.We 
expect that the finite-temperature superconducting tran- 
sition will be of the conventional KT type, triggered by 
the unbinding of integer vortices. 

It turns out that the triangular lattice of coupled grains 
is special in its high degree of order at low temperature. 
The square lattice for d ± id grains is at the opposite 
extreme, as shown in the following section: there is no 
coupling between the Ising variables and the supercon- 
ducting phase. The zero-temperature degeneracy for the 
square lattice is thus given by a global U(l) for the su- 
perconducting phase times a factor of 2 per site from the 
Ising variables. 

The case of p ± ip superconductors on the square lat- 
tice is quite complex and is shown below to be connected 
to the existence of a Z 2 topological order. We find that 
p±ip superconductors on the square lattice support lo- 
calized 7r-vortices with half the enclosed flux of an or- 
dinary superconducting vortex, and these vortices drive 
the superconducting transition. 
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FIG. 3: A typical ground state configuration of Ising variables 
(open or filled circles) for p + ip superconductors. Note that 
fixing one row and one column determines the rest of the 
lattice via the constraint of 0, 2, 4 up-spins per square. 



V. SQUARE LATTICE: DECONFINEMENT 

Square arrays of both p+ip and d+id superconductors 
show dramatically different behavior from the triangular 
lattice discussed above: both have no order of the Ising T 
variables even at zero temperature when only the Joseph- 
son coupling is considered, i.e., only single-Cooper-pair 
processes are included in the model. 

The gauge flux $ a that pierces an elementary plaquette 
in the square lattice is given by 



k 

$a = + (72 + 03 + (T 4 ) 



xi 1. 



(15) 



Again k = 1 for p + ip and k = 2 for d + id. For d + id the 
flux generated by the Ising variables is | Ylo CTi > which is 
always an integer. Since integer gauge flux can always be 
gauged away, the problem decouples into a square-lattice 
XY model and a non-interacting Ising model. 

The case oip+ip arrays is more interesting. In this case 
$ a = i Gi and hence can be fractional. Indeed, for 0, 
2, or 4 up spins per square $ a is an integer, while for 1 or 
3 up spins $> a is \ plus an integer. In the ground state, 
the Ising spin configurations must satisfy the constraint 
that there are even number of up spins in every plaquette. 
In a system with open boundary conditions and size M by 
N, the number of such Ising configurations is 2 M x 2 . 
The first factor corresponds to the number of ways to 
fix all the spins on one horizontal line, and the second 
factor the number of ways to fix all the spins on one 
vertical line. Once those spins are fixed, the remaining 
spins are uniquely determined by the constraint (Fig.Pjl. 
This implies that, although the above system appears 
two-dimensional, its ground-state entropy grows as the 
system's linear size rather than its area. 

The ground-state constraint of 0, 2, or 4 up-spins per 
square corresponds to a degenerate case of the eight- 
vertex model^, which can be mapped onto the six- vertex 
models using a duality transformation^. The point 
described by \ Y^n °i integral is quite special: it di- 
vides a region of the phase diagram with extensive (two- 
dimensional) entropy from a region with finite (zero- 
dimensional) entropy. 



The ground states for p+ip superconductors on the 
square lattice are the same as those of the Ising lattice 
gauge model 



H = -K 



i,3 



(16) 



since the product around a plaquette takes the value +1 
for 0, 2, 4 up spins, and —1 otherwise. A plaquette with 1 
or 3 up spins generates a gauge flux. Whether such defec- 
tive plaquettes are entropically confined is equivalent to 
whether the the 7r-gauge flux is confined in the gauge the- 
ory when K — > 0. Our previous degeneracy counting im- 
plies that the 7r gauge flux would completely deconfincd 
if the entropic term were the only contribution to its free 
energy. In reality, in addition to the entropic contribu- 
tion, there is a phase stiffness contribution to the inter- 
action between two fractional gauge flux. Indeed, when 
<f> a is a fraction, it cannot be screened by integer vortices 
n a . As the result, two fractional fluxes <E> separated by a 
large distance R cost free energy ~ (nK/4) log-R, where 
K is the renormalized phase stiffness. This logarithmic 
interaction binds half-integer vortices at low tempera- 
tures. 

The half-integer vortices are relevant at the ordinary 
(integer vortex) XY transition because their scaling di- 
mension is four times smaller than that of the integer 
vortices, which are marginal at the transition. It fol- 
lows that as the temperature is raised above a critical 
value, the half-integer vortices unbind and the supercon- 
ductivity is lost. The superconducting phase below this 
temperature has "nematic" order in that there is long- 
range order in the square ip 2 of the Cooper pair operator, 
rather than in ip. 

A defective plaquette can be viewed as a 7r-flux bound 
to a fractional Ising charge. The Ising charge associated 
with a given plaquette is defined as 



Qi = 



(17) 



where □ is the number of up spins in the plaquette. 
Any plaquette satisfying the 0, 2, 4 constraint has an 
integral Qj. An isolated plaquette with 1 or 3 up spins 
has half- integer Qj. Since the Josephson energy is only 
sensitive to whether the Qi is integer or half-integer, 
one should regard Qi as defined modulo integer since 
two half-integer Qj cancel. Local constraints like (|15|) 
were used in^l to generate a gauge theory, but required 
a somewhat artificial charging energy that is associated 
with plaquettes, not with single grains or with long-range 
interactions. We have shown that geometric phases in 
Josephson junction arrays naturally generate such a pla- 
quette constraint, and in turn a classical two-dimensional 
Z 2 gauge theory. 
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VI. HONEYCOMB LATTICE: CRITIC ALITY 

We now study coupled p ± ip and d ± id grains on 
the honeycomb lattice and show that the essential low- 
temperature physics is considerably modified from the 
triangular and square cases studied in the preceding sec- 
tions. This model will be shown to have an extensive en- 
tropy down to zero temperature with no macroscopic or- 
der of the Ising sign variable. Hence even though locally 
(i.e. in each superconducting grain) T-breaking super- 
conductivity is favored, globally there is no time-reversal 
symmetry breaking; there is instead a critical state at 
zero temperature. 

For the honeycomb lattice, the gauge flux in each 
hexagon is given by 

k 

$a = 77(0-1 + 0"2 + 03 + (7 4 + O5 + (T 6 ), mod 1. (18) 
6 

Hence a plaquette with zero, three, or six up-spins has 
integer <E> a , and does not cause any frustration in the 
Josephson energy. Thus the lowest energy Ising spin 
configurations satisfy the 0, 3, 6 constraint in every pla- 
quette. The number of such configurations can be found 
exactly in the thermodynamic limit by a mapping 28 (sec- 
tion IVIII) onto the three-color model solved exactly by 
Baxter—. The result gives an entropy per hexagon equal 
to 0.379114fcB- The critical properties of this three-color 
model are briefly reviewed in section VII. 

Ising configurations that violate the 0, 3, 6 constraint 
lead to fractional gauge fluxes. The number of Ising con- 
figurations consistent with a distribution of fractional 
gauge fluxes is different from that of the ground state 
degeneracy discussed above. The logarithm of the ratio 
between these degeneracies gives rise to a entropic contri- 
bution to the interaction free energy between fractional 
gauge fluxes. To compute such interaction we need to 
consider a correlation function that does not have an ana- 
logue in the three-color model: this is the main subject 
of section VII. Here we explore the physical consequences 
of the entropic interaction. 

To be more precise let us return to Eq. (|12f> . For a 
fixed distribution {$ a } (-1/2 < $ a < 1/2) the entropic 
contribution to the dimensionless free energy is given by 



-AF ({*„}) _ ^{°*> H« H $a S a i 



(19) 



Here the function 8(x) is 1 if x € Z, otherwise. After 
the summation over the Ising variables, Eq. I|12|) becomes 



Z cx£ 



-S({n a ,$ a }) 



(20) 



where 



S = A c ^(n a - $ a ) 2 +J2( n a- ®a)G ab (n b - $ 6 ) 

a a,fc 
+Af({*a}). (21) 



The question is what effect, if any, does AF({$ a }) have 
on the thermodynamics of the vortices? 

Section IVIII shows that the probability e~ AF associ- 
ated with two $q = ±1/3 defect plaquettes of opposite 
fractional flux falls off with distance approximately as 
r -4/9_ This implies a free energy cost of 



AF = 4 $ a log I 



(22) 



Imagine the stiffness K is tuned so that the system is 
exactly at the critical point where integer vortices pro- 
liferate if fractional vortices are strictly forbidden. Then 
the vortex fugacity operator is marginal (has scaling di- 
mension 2). From the Gaussian spin- wave theory vortices 
of fractional charge 1/3 have scaling dimension A = 2/9. 
However, this analysis ignores the entropic contribution 
to the free energy. Including that effect from Eq. im- 
plies the total scaling dimension of fractional 1 /3 vortices 
at the (integer vortex unbinding) KT point 



A = 2/9 + 2/9 = 4/9. 



(23) 



This indicates that fractional vortices are in fact rele- 
vant, and most likely it is their unbinding that drives 
the superconductor-normal XY-type transition. The su- 
perconducting phase is then a "triatic" with long-range 
order in the cube of the Cooper pair operator ip 3 rather 
than ip. 

There are several corrections in realistic systems 
(quenched disorder, multiple-Cooper-pair processes, etc.) 
that will split the extensive zero-temperature entropy, as 
discussed briefly in the final section and in a forthcoming 
paper by Castelnovo, Pujol, and Chamon— However, 
these effects can be expected to onset at temperatures 
much less than the Josephson coupling (kT <C Ej), so 
experiments on the honeycomb lattice at the lowest ac- 
cessible temperatures are likely to see an extensive en- 
tropy and critical correlations of the time-reversal order 
parameter. The following section obtains these critical 
properties using the connection to the three-color model. 



THE T = CRITICAL THEORY ON THE 
HONEYCOMB LATTICE 



VII. 



As shown in section IVII the ground state configura- 
tions of the Ising spins have 0, 3, or 6 up-spins on each 
hexagon. The entropy per site of this "0, 3, 6- model" was 
earlier shown^^S to be equal to that of the "three-color 
problem" solved by Baxter. In the three-color problem, 
one assigns one of three colors (say red, green, blue) to 
each bond of the honeycomb lattice, with the require- 
ment that all three colors meet at each site. A quantity 
called chirality can be defined in the following way. At 
each site there are three bonds each being colored differ- 
ently. Start with the red bond and turn clockwise: if the 
sequence of colors encountered is RGB we assign chirality 
cr = +1 to that site. If the sequence is RBG we assign 
a = — 1. In this way each coloring pattern is uniquely 



6 



associated with an Ising configuration. The 0, 3, 6 con- 
straint ensures that one returns to the original color after 
going around a plaquette (hexagon) . Since the colors are 
uniquely determined from the Ising variables once one 
bond is fixed, the entropy of the two models differs by 
a global factor of 3. Baxter has shown that the entropy 
per hexagon is 1 - 



Width 



Estimated S/o 



X, 



log 



1-34-67-9 



0.37905. 



(24) 



The three-color model is an example of zero tempera- 
ture models that have critical correlations^, and many 
of its critical properties are known exactly. To determine 
these properties it is useful to represent the three-color 
model as a two-component height model (on the dual, tri- 
angular, lattice) as follow a 1 ^' 1 ^? 2 . Consider a given site 
on one sublattice (say A), and the three colored bonds 
stemming from it. If we cross a given bond in a clockwise 
fashion with respect to the site, the following increment 
to the high variable is assigned: 



Ah = (1,0) red 

= (-1/2, V3/2) blue 
= (-1/2, -V3/2) green 



(25) 



This assignment is self-consistent since the three-color 
constraint ensures the total height increment is after 
one completes a loop. The height increment gets a factor 
of —1 if the crossing is counterclockwise with respect to 
sublattice A. The long- wavelength effective theory for 
this problem is given by 



F = 



^(Vh) 2 + y(h)) ,/-,-. ,.>(i) 



where the potential V has the periodicity of a honeycomb 
lattice of side 1 / y/3. The stiffness K is determined by the 
requirement that the locking potential V be marginal. 

In this theory the allowed vertex operators have the 
form O a (x) = exp(iaihi(x)) where a lies on the recip- 
rocal lattice of V. The scaling dimension A of such 
an operator, which determines its correlations through 
(6 a (0)6 a (x)) ~ r- 2A , is given by A = a 2 /(3ir). Know- 
ing the critical theory of the three color problem it is pos- 
sible to compute the entropy associated with certain lo- 
calized defects in the color assignment. These defects ap- 
pear as vortex operators; among them the most relevant 
one has scaling dimension A = 1/2. This information 
allows us to determine the full critical three-color theory. 
This critical theory l|2t)|l has a hidden SU (3) symmetry^. 
The field-theory representation of the SU (3) symmetry 
of i|26|) is given in3£. 

Unfortunately knowledge about the three-color model 
does not enable us to compute the entropic contribution 
to the interaction between fractional gauge fluxes as in 
Eq. I|22|) . This is because a single plaquette with the 
wrong number of up Ising spins spoils the mapping to 



Open 



Periodic 



1.05 
0.703 
0.590 
0.534 
0.502 
0.480 
0.465 
0.454 



0.462 0.210 



0.403 0.231 1.90 



0.390 0.235 1.91 



Theory (oo) 0.379114 0.379114 2 

TABLE I: Results of numerical transfer-matrix calculations: 
estimates of T = entropy and correlations. The correlation 
between two defect plaquettes of opposite sign (e.g., one with 
2 spins up and one with 4) falls off as r~ 2Ad , and the chirality 
correlation falls off as r~ 2Ac . 



the three-color model. For example, on circling a de- 
fect plaquette with only 2 Ising spins, the colors wind up 
cyclically permuted. To find the entropic interaction be- 
tween a pair of fractional vortices with opposite vorticity, 
we numerically computed the ratio in Eq. (|19fl . 

We have performed numerical transfer-matrix calcu- 
lations on strips and cylinders of up to 9 hexagons 
width/circumference. The entropy for the ground state is 
calculated as follows: configurations allowed by the con- 
straint (0, 3, or 6 up-spins on each hexagon) are assigned 
(3E = 0, while forbidden configurations have f3E = oo. 
Let the length of the strip be N: the entropy is given by 



S = \ogZ = \og(T N ) = Nlog\ 1 



(27) 



where Ai is the largest eigenvalue (assumed nondegener- 
ate) of the transfer matrix T. The entropy per hexagon 
of the 2D system, if extensive, can then be estimated as 



S/o = lim 



logAipy) 

w ' 



(28) 



where W is the width of the strip. Table I shows the re- 
sulting estimates with periodic and open boundary con- 
ditions for strips of various widths. 

The transfer matrix can also be used to calculate the 
entropic interaction between a pair of plaquettes with op- 
posite flux on a cylinder. These defect plaquettes must 
be on the same sublattice for a nonzero correlation. We 
use the usual procedure^ for extracting scaling dimen- 
sions at a critical point from calculations on the cylinder, 
based on the hypothesis of conformal invariance. The 
relative probability to have positive and negative defect 
plaquettes at fixed locations, described by transfer ma- 
trices T + and T_ , separated by n defect-free rings around 
the cylinder, is 



p(n) 



Z(n) 



lim 

N^oc 



Tr ' Fj r T n T 



(29) 
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If this function in the continuum limit is a correlation 
function of a critical theory in 2D, then conformal in- 
variance predicts that correlations along the cylinder at 
large distances r ^S> L, L the circumference of the cylin- 
der, should fall off as 

p(r)~e-^, C=2^. (30) 

As seen in Table I, this characteristic dependence is in- 
deed observed for both the fractional-vortex-defect cor- 
relator and the Ising correlator. 
The numerical results suggest that 

-AF(r a ,r b ) 1 /oi\ 

\r a -n\^- (31) 

This implies that the fractional vortices corresponds to 
an operator with scaling dimension A w 0.235 in the 0, 
3, 6 model. Since the lowest two scaling dimensions that 
appear in the three-color model are A = 1/2 and A = 2/3 
which are very different from 0.235, it is apparent that 
the critical theory of the 0, 3, 6 model has a different 
operator content than the three-color model. 

Such a critical theory must have the same correlators of 
height variables as the three-color model, and yet have 
a different set of defect operators, including the defec- 
tive plaquette with A w 0.23. One possibility is the 
"orbifold" dual of the above three-color model theory ob- 
tained by identifying heights that differ by a Z3 transfor- 
mation h\ + ih% — * exp(2Mrn./3)(/ii + ih^)- The physical 
significance of this identification is that the fractional 
vortices, around which the height lattice is rotated by 
27r/3, are now allowed defects. On the torus, there are 
now different sectors of the theory corresponding to how 
the color scheme is shifted upon passing from one bound- 
ary of the torus to another. We find that the orbifold 
"twist" operators that connect these different sectors 3 ^ 
have scaling dimension 2/9 and may correspond to the 
fundamental plaquette defects in the 0, 3, 6-model. How- 
ever, this identification is quite indirect at present and 
should be regarded as only one of several possibilities 
for the critical theory. This orbifolded SU(3) theory 
with c = 2 has previously been considered as a hypo- 
thetical model of a quantum Hall edge 3 ^, since the non- 
orbifolded SU(3) theory describes several real quantum 
Hall edges2§i2£ in the chiral Luttinger liquid approach 38 . 

Finally, Eq. (|31|l indicates that fractional vortices have 
already become unbound when the integer vortex pairs 
disassociate. To see that we imagine forbidding the frac- 
tional vortices and tuning the stiffness constant K so 
that the integer vortex just about to unbind. At that 
point the vortex fugacity operator is marginal (has scal- 
ing dimension 2), and 1/3 vortices have scaling dimension 
Aj = 2/9. Adding the entropic interaction (Eq. i|31|)1 this 
yields the total scaling of the fractional vortex fugacity 
operator 

A = A ; + A e = 2/9 + 2/9 = 4/9. (32) 



This indicates that fractional vortices are in fact rele- 
vant, and most likely it is their unbinding that drives the 
superconductor-normal XY-type transition. The super- 
conducting state has "triatic" order in ip 3 rather than in 
tp, as in the ground state of the antiferromagnetic XY 
model on the kagome lattice 18 . 



VIII. SUMMARY AND OTHER DIRECTIONS: 
DYNAMICS AND QUENCHED DISORDER 

Our main goal in this paper was to demonstrate that 
due to the coupling of the order parameter and the under- 
lying crystalline lattice, Josephson junction arrays made 
of time-reversal-symmetry- breaking superconductors ex- 
hibit very much richer physics than their s or d coun- 
terparts. Phase-sensitive experiments such as those in 
Refi&i should be able to detect the variable tunneling 
phases caused by a T-breaking order parameter. Regu- 
lar arrays of such grains give rise to realizations of Higgs- 
gaugc models, and are new examples of topological order 
with deep connections to frustrated magnets. The order- 
ing of the time-reversal order parameter and the onset of 
superconductivity depends intricately on the geometry of 
the arrays. 

One can expect, on a qualitative level, that the ex- 
istence of the nontrivial phases studied above will also 
modify the behavior in random polycrystalline samples. 
The behavior in the random case is likely quite com- 
plex since random Josephson coupling, even without the 
gauge phases we have studied, leads to the "dirty boson" 
problem 3 ^^. Allowing random local orientation of the 
order parameter for d superconductors has been argued 
to induce a superparamagnetism effect^ that will coexist 
with the spontaneous orbital currents in the d+id case. A 
more tractable model would be to allow a quenched time- 
reversal-symmetry breaking field coupling to the Ising 
variables, but keep uniform the strength of the Joseph- 
son coupling. 

The possibility of a direct coupling between the Ising 
variables at higher order in the Josephson tunneling is 
studied in a forthcoming paper—. This model has phase 
transitions from the critical phase with no additional in- 
teraction to ordered ferromagnetic and antiferromagnetic 
phases. Its dynamics are of special interest, as the three- 
color model has previously been studied as a model sys- 
tem for glassy dynamics*. Our realization of the three- 
color model in arrays of complex superconducting grains 
may prove to be experimentally less challenging than 
recent efforts^ 3 to build frustrated XY models on the 
kagome lattice using s superconductors and a frustrat- 
ing magnetic field. 
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